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Abstract 

An initial-boundary value problem for the ID self-adjoint parabolic equation on the 
half- axis is solved. We study a broad family of two-level finite-difference schemes with 
two parameters related to averagings both in time and space. Stability in two norms 
is proved by the energy method. Also discrete transparent boundary conditions are 
rigorously derived for schemes by applying the method of reproducing functions. Results 
of numerical experiments are included as well. 
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1 Introduction 

In many applications, a problem of solving partial differential equations 
in unbounded domains arises. A number of approaches to the problem is 
developed mainly associated with the statement of additional boundary con- 
ditions on artificial boundaries [l]-[6]. The conditions are called the (exact) 
artificial/non-reflecting/transparent boundary conditions provided that they 



are satisfied by the solutions of the original problems in unbounded domains. 
For definiteness, we exploit the last name (TBCs). For parabolic evolution 
equations or for the Schroodinger equation, the TBCs are integro-differential 
relations along the artificial boundaries. Their adequate discretization is non- 
trivial since it can produce significant reflections from the boundaries and even 
instability in computations as well as create difficulties in rigorous proofs of 
stability of the resulting numerical method. 

An alternative approach suggests to implement the idea of the TBCs on 
the mesh level. Namely, first one consider a discretization of the problem on 
an infinite mesh in the unbounded domain (which is not practical because 
of the infinite number of unknowns). Its solution is restricted to the finite 
mesh by deriving a mesh counterpart of the TBC on the artificial boundary. 
One version of this approach is associated with the derivation of discrete 
TBCs requiring to solve analytically model mesh problems on infinite grids. 
This approach had worked well by the complete absence of reflections from 
artificial boundaries and reliable stability of computations in practice as well 
as by clarity of the mathematical background and rigorous proofs of stability 
of the resulting mesh method in theory. Such an approach is developed in 
detail for ID time-dependent Schrodinger equation in [7] -[11] and also used 
for ID parabolic equations in [1, 12, 8, 6, 13]. 

In this paper, the approach is developed for the ID self-adjoint parabolic 
equation on the half-axis. We study a broad family of two-level finite- 
difference cr-schemes with averaging in space with a weight 6. We prove 
its stability by the energy method and rigorously derive the discrete TBC by 
applying the method of reproducing functions. Notice that both points are 
not so well developed in some other papers. Similarly to [11], this allows one 
to cover in a unified manner a collection of particular schemes: the standard 
scheme without averaging (9 = 0), the linear finite-element method (0 = |), 
scheme of higher order of accuracy (for constant coefficients, 6 = j^) and a 
vector scheme on a four-point stencil (6 = j). The results generalize those 
obtained for 6 = in [13] but here we prove stability in two (not one) energy 
norms, the rigorous derivation of the discrete TBC is notably different and 
results of numerical experiments are included as well. 



Notice that the results for the family of finite-difference schemes can be 
enlarged for the 2D (or multi-D) case exploiting the technique from [14, 15]. 



2 An initial-boundary value problem and a finite- 
difference cr-scheme with averaging in space and an 
approximate TBC 

We consider the one-dimensional parabolic equation 

P^ + Au = f, An -- (b-j + cu (1) 

for x > and t > 0. Its coefficients satisfy the conditions p(x) ^ p > 0, 
b(x) ^ v > and c(x) ^ for x > 0. 

Equation (1) is supplemented with the following boundary condition, the 
condition at infinity and the initial condition: 

u \x=o = 9{t)i u(x,t) — > as x — > oo, for all t > 0, (2) 
u\ t=Q = u°(x) for x > 0. (3) 

We assume that the coefficients become constants as well as / and u° vanish 
for sufficiently large x ^ Xq (for some Xq > 0) 

p{x) = poo > 0, b(x) = boo > 0, c(x) = Coo > 0, f(x,t) = 0, u°(x) = 0. (4) 

An integro-differential TBC satisfied by the solution to this problem can 
be written in the Dirichlet-to-Neumann form 

dx K ' \l boo V^dtJ v J ^t^e v ; 

for t > and for any X ^ Xo; other equivalent forms are also known. The 
TBC is nonlocal in time; recall that the involved operator 



^1/2 / s 1 d ... d9 

defines the classical left-hand Riemann-Liouville time derivative of order | on 
the half-axis [0, oo). But we do not exploit this TBC explicitly below. 



We fix some X > X , set Q = (0, X) and introduce a nonuniform mesh 
Uh,oo m x on [0, oo) with nodes = xo< ••• < xj = X < . . . and steps 
hj := Xj — Xj-i such that hj ^ X — Xq and hj = h = hj for j ^ J. We set 



Xj-i/2 := Xj and h J+1 / 2 ■= j+ 2 j+1 • Let u) hj00 := w/ l , 00 \ {0}, co h := {xj}j =Q 
and Uh '■= {xjYjZi- Let Ho(uJh) be the space of functions on mesh uJh that 
equal for xq = 0. 

We define the backward, modified forward and central difference quotients 
in x 



Ai, Aii -Li /9 ^Aii 



ij ' 4 i+l/2 z/ ij+l/2 

together with averaging operators in £ 



— ^ + 

+1/2 

kn ... , , „ hn- 



s x W 3 := ^-Wj + 

zn j+l/2 z > l j+l/2 



Se Wj := e-r^—Wj-i + (1 - 2(9) jy ? + 

hj+l/2 ' %+l/2 

and a more general mesh counterpart of multiplication by a mesh function 
CffMWi := + (1 - 20) (s^)^- + ^ir ±L ^ J+1 Wi +1 , 

hj+l/2 "i+1/2 

see [11]; clearly C#[l] = sg- Notice that sqWj = s d Wj + sjWj, where 

s Wj := 0Wj_i + (\ - d)Wj, s+Wj := (± - 0)PFj + 0W>+i. 

We also introduce a nonuniform mesh in t on [0, oo) with the nodes = to < 
■ ■ ■ < t m < . . . such that t m — > oo for m — > oo and the steps r m := t m — t m -\. 
Let uj t := uJ T \ {0}, uf M := {t m }m=i- We a ^ so define the backward difference 
quotient, the mean value with the weight o (independent of the meshes) and 
the backward shift in t 

d t $ m ■= , $W m : = cr$ m + (1 - cr)<I> m_1 , <l m := $ m_1 . 

We study the following finite-difference scheme, which is weighted in t and 
averaged in x, on a finite mesh with an abstract approximate TBC for the 



initial-boundary value problem (l)-(4) 

Ce[p h ]d t U m + A h U^ m = F m on u h x W T , (6) 

= g(t m ) for m > 1, (7) 

bood x U^ )m + hjSe( Poo dtU + cboOW)? = 6oo<S w U7 for m ^ 1, (8) 

C/° = U° h on 57 h (9) 

with the operator AhW := — d x (bhd x W) + Co[ch]W and the functions 

P/y = p(Xj-i/ 2 ), b hj = b{Xj-i/ 2 ), C hj = c(x J _i/ 2 ), Ff = f(Xj,t m ) 

and Ujlj = u°(xj) (for simplicity, for continuous p, 6, c, / and ito). Thus 
Uflj = 0; we also assume that U^ = 0. Here S m is any linear operator acting in 
the space of functions given on the mesh w^U{0}, and := {U% . . . ,Uj). 

Now we discuss the approximate TBC, i.e., the boundary condition (8). 
Let an equation 

eh 2 



d 



S m UJ for m ^ 1 (10) 



serve as an (abstract) approximate TBC for (5) at the node xj\ here we 
have discretized |^ with weight in t and symmetrically in x. We first write 
down equation (6) on the mesh Uh U {xj} and apply it at the node xj only 
in order to eliminate the values Uj\ involved in the left-hand side of (10). 

o ^ 

Namely, since d x Wj = d x Wj + kd x d x Wj, taking into account (4) we get the 
boundary condition (8). Importantly, the boundary condition (8) for S = is 
the natural approximation of the Neumann boundary condition for this finite- 
difference scheme for x = X. Such an approach was implemented earlier in 
[9, 13, 11]; it reliably leads to the computationally stable form of discrete 
TBCs (in contrast to some other approaches). 

The corresponding three-point system of mesh equations for a vector 
{C/^j^o of the solution values on the upper level, has the form: 



for 1 < j ^ J - 1, 

cc = g(t m ), o^jUf^ + (2^_ lfJ - fi m0 )uj 

= «?-l,J^l 1 + (2^_ ljJ -M(m-l)l)^ l - 1 + E Mi(m-0^. (11) 

0<Z<m-2 

compare with [11]. Here the coefficients are given by formulas 

- 61 — 1" <™jC fc j I - cr— , (j aJ -(--V)l — h ahjChj I + a— 

for 1 ^ j ^ J; in particular, for j = J, these expressions become more simple 

m n ( hPoo , \ boo nm /l n\ ( ^Pco , \ ^00 



Equation (11) is written assuming that the operator in the approximate TBC 



has the form S m UJ = £ Vl(m-l)Uj. 

i=i 

3 Stability of the finite-difference scheme on the finite 
mesh with the approximate TBC 

We consider the stability problem for the finite-difference scheme (6)- (9) 
with respect to the initial data U®, the free term F and a perturbation in the 
boundary condition (8) and take g(t) = 0. We need to introduce several mesh 
counterparts of the L 2 (Q)-inner products 

j-i j 

M = E VWW, (Y, w)~ h = E VjWjhj, 

3=1 3=1 

(Y, w)~ h = (V, w) Uh + VjwA 



II ' llw h > 


' 1 1 0J h 5 


' \\u> h 









and the corresponding norms 
given on u>h or belonging to Ho(uJh)). 
We introduce a bilinear form 

(U, W) Ce[ „ h] := (ColxhP, W) Uh + j^j/ij for U, W € H Q (uJ h ). 

According to [11], it is symmetric for ^ \ and generates a norm HWH^r^i = 

1 /2 

(W, W)^,^ for Xft > (or a seminorm for x/j ^ 0). Moreover, an inequality 



V^IWk < ll^lb eW < V / (l + 4max{-^,0})p||W|k (12) 



holds for any W G H^uju) and ^ |, where c,? = 1 — 4max{#,0} and 
p = max p hj . 

We also introduce mesh counterparts of the norms in L 2 (0, tu) and L 2 (J1 x 
(0,t M )): 

M 



l*lk : = (E(* m ) 2 ^) 



m=l 



and also set 



\C 9 [K h ] 



M 



x 1/2 

7 ' 








WujhWujIj 






\cjr M for ^ 0. 







Proposition 1 Let U be a solution to the finite- difference scheme (6)- (9) 
with a generalized boundary condition (8): 

bood x uf )m + hjs e (p h d t U + c h U^yj = &oo<S m ir; + G m for m > 1, (13) 



where G is given on uf . Let the operator S satisfy an inequality 

M 

(S m $ m ) <5> {a)m r m ^ for any M > 1 



(14) 



m=l 



/or an?/ function $ gwen on a; T swca £/m£ $° = 0, where $ m = {$\ . . . , $ m }. 
Then, for a ^ | and < \, the first energy bound 



m ax{ m^ M ll^ m ||c 9 H,V2||C/||(i)} ^ II^ILm 



M 



'cop 

— m=l 



^ || F (0)i 



\uj h T m 



+ ^-(l|F (1 »ll sV + Vx||G|k) 



(15) 



ao/ds /or an?/ M ^ 1 and /or an?/ decomposition F = + d x F^ such that 
F^\j=j = 0, with K a := 2(o~ + |1 — i/ere £ae norm ||?7||m «s swca £aa/: 



lic/||?i, 



" - m^d t u\\ 2 c, lnK , + llv^^'ll^ + l|t/ ( 



cr) || 2 



Tae bound holds also in the case 6 = | provided that F^ = ( one has to 
drop the summand with F^>). 

Consequently, for a ^ \ and 9 ^ \, the scheme has a unique solution. 



Proof. We take the 



1 JOJh 



-inner product of equation (6) and a function 
W G Ho(uJh), sum the result by parts (using the second assumption (4)) and 



obtain 

{ColptAdt^WV + (b h d x U^ m : d x w)„ + (Ce[c h }U^ m : W 



= b 00 (B x U^ m WT + (F m ,W) u}h . (16) 

Choosing W = U^ m and applying the boundary condition (13) and other 
assumptions (4), we get 

(d t U m , U^ m )c e[Ph] + (b h d x U^ m ,d x U^ m )z h + (U^ m , U^ m ) Ce[ch] 
-b^ (S m VJ) Uf )m = (F m , U {a)m )u h + G m uf )m for m > 1. 

We multiply the result by r m and sum up it over m = 1, . . . , M. Applying 
the formula = £/ (1/2) + (cr — h)rdtU, we obtain the first energy equality 



M 

+\\v^\? CAChYM - b„ y: (*s m u-;) u? )m T m = m\ 2 c M + i {i)m (") 

771=1 

for M ^ 1, where 

j(l)M . = ^ I (^m ^Wm)^ + Qmjji 



M 

r((j) m 

' J 

m=l 

For F = F^ + d x F^\ we sum the result by parts and derive the following 
bound 

M 

^ J- \\F^ m \U (MIlCHk + |1 - alH^IL) r m 



771=1 



+(11^113^ + ^l|Glk)lia*f("|| a!|r 

M 

< (M + II - *l) E ll^ ( ° )m |l^ max ||[/ 



m=l 



Using the left-hand inequality (12), conditions a ^ \ and (14) and applying 
the standard argument lead from the first energy equality (17) to bound (15). 
It is well-known that such a bound implies the existence and uniqueness of a 
solution to the finite-different scheme. 



Remark 1 In the case 9 = \, one can generalize bound (15) and next 
stability bounds for ^ as well, see [11]. 

We also define a symmetric bilinear form 

c UM (u : w) = (b h d x u : d x w)z h + (c e [ Ch \u, w) Uh + c^u^Wjhj 

1 /2 

for U,W G i7o(cJh) and derive stability in the norm ||VF||£ = £J hg (W, W). 

Proposition 2 Let U be a solution to the finite- difference scheme (6)- (9) 
with the generalized boundary condition (13) instead of (8). Let the operator 
S satisfy an inequality 

M 

(S m $ m ) d t <5> m T m ^ for any AO 1 (18) 

m=l 

for any function $ given on uJ T such that $° = 0. Then, for a ^ -| and 9 < \, 
the second energy bound 

max{ max 11^11^,^11/711^} < \\U° h \\ c + J— \\F®\Ur 

. M M 

V (1) °k + E P^ (1) kTm + >/*(|G°| + E l^> m )] (19) 



m=l m=l 



holds for any M ^ 1 and am/ decomposition F = F^ + d x F^ with F^\j=j = 
0. ii/ere 

M 



TO=1 



rm\\2 
\C e [p h ] 



Tn 



The bound holds also in the case 9 = | provided that F^ = 0. 
Consequently, for a ^ | and ^ 7, £ae scheme has a unique solution. 

Proof. We choose = c^£7 m in (16), apply the boundary condition (13) 
and assumptions (4) and get 

(d t u m , d t u m ) Ce[ph] + c UM {u& m , d t u m ) 

-600 (S m V n j l )d t Uy = {F m ,d t Uj) ujh + G m d t UJ for m > 1. 



We multiply the equality by r m and sum up it over m = 1, . . . , M. Applying 
again the formula = U^ 1 ^ + (a - \)rd t U ', we obtain £/ie second energy 
equality 

M M 



m=l 



\ l T + 1||[/M||2 



M 



7i« 



777 = 1 



-600 V (<S™U?) ^c/yv TO = i ||t/, ||^ + / (2) 



m 1 1 2 _2 

M 



(20) 



771=1 



for M ^ 1, where 



M 



/( 2)M . = £ [ (F m 5 + gm^] ^ 



m=l 



For F = FW + dxF®, we sum the result by parts in t and x and get 



7 (2)m = ^( J p(o)m j d^"% ft T m - (F® m ,d x ir 



777=1 



TO=M 
777=0 



M 



M 



<||F( )|U T ||9 t C/|L- + 2[||F« | 



777=1 



777=1 



M 

M 



M A/ 

+ £ ll^ (1) lkr m + + £ \d t G m \r m ) 



777=1 



777=1 



max ||5 x 7/ m || s?i . 

0<m<M 



Using the left-hand inequality (12), conditions a ^ \ and (18) and applying 
the standard argument lead from the second energy equality (20) to bound 
(19). 



4 Stability of the finite-difference scheme on an infinite 
mesh 

In order to construct and study the discrete TBC, we first turn to the 
finite-difference scheme on an infinite mesh for the original problem (l)-(3) 



on the half-axis 

Ce[p h ] d t U + A h U^ = F on oj h}00 x uT, (21) 

U^ = g(t m ) formal, (22) 

U° = U° h on uJ Koo . (23) 

Assumptions (4) are supposed to be fulfilled. Let g(t) = and U^\j=o = 0. 

We introduce the Hilbert spaces Hh and Hh (mesh counterparts of L 2 (R + )) 
consisting of functions W given on the meshes respectively uJh,oo (and with 

oo 

Wo = 0) and U0h,oo and such that ||W||f = Wf < oo, equipped with the 

3=1 

inner products 

OO oo 

(V, W) Hh := Y, VjWjh j+1/2 , (V, W) Sh := £ V 3 W 3 h 3 . 

3=1 3=1 

Since hj = h for j ^ J, the conditions ||VK|| 2 9 < oo, ||W||jy < oo and 
||W||~ < oo are equivalent. 

H h 

We define symmetric bilinear forms [11] 

(V : W) Ce[ « hi00 := (CeW]V : W) Hh: 
£hJV,W) = (b h d x V,d x W) 5h + (Co[c h ]V,W) Hh 

for V,W G H h . They generate norms ||V^|| Cfl [« h ], o = W w )g(k]^ for x ^ = 

l li 

Ph (a seminorm for Hh = Ch) and ||W||£ ff = jC^ hg (W,W). Moreover, an 
inequality 

y/cop\\W\\ Hh < 11^11^^,00 < V / (l + 4max{-^, 0})p||^|| Hfc 

holds for all W G # ft and 6 < ± [11]. 

We also define the mesh counterparts of the norm in L 2 (K + x (0,t/y/)) 

II • \\ H ^ '■= II II ' lUJUlf II • llj^ T := II II • ll#Jkv 

Proposition 3 Let F = F<® + d x F^ with F^ m G H h and F^ m G H h 

for any m ^ 1 and G ii/^. Then, for a ^ an<i# ^ j, there exists a unique 
solution U m G Hh, for all m ^ ; to the finite- difference scheme (21)- (23), 



and, for a ^ | and < \, the first energy bound 



^ M 

^ WuftWcwfio + — = E ll^lk^ + \/^ H F(1) b- (24) 



/cep ^— ' V ^ 

— m=l 



/io/ds /or an?/ M ^ 1 . Here 

2 



U 



|2 
l(l),oo 



v /(o--i)r||^C/|| CeW)00 +1111^11^ 



|2 



T/ie bound holds also in the case 9 = | provided that = 0. 

Proof. We extend ^ and C^p/J up to operators acting in by setting 
(AhW)o := and (C^p/jV^o := 0. By virtue of assumptions (4) and the 
property hj = h for j ^ J, the operator ^ is bounded Hh- Moreover, 
Ah = A* h > since 

(AhW,V) Hh = C Hm {V,W) (25) 

for any W, V G i^. To establish equality (25), one can first transform the 

h 

finite sum {AhW)jVjhj + i/2 by summing by parts (compare with the deriva- 

tion of equality (16)) and then pass to the limit as j\ — > oo using the property 
lim^oo Wj = for W € H h (as in [13] for = 0). 

Now we rewrite equation (21), together with the homogeneous boundary 
condition (22), as an operator equation in Hh 

CeiPh] dtU + A h U^ = F on u T . (26) 

Thus 

(Celp h ] + ar m A h )U m = (C e [ Ph \ - (1 - a)r m A h )U m + r m F m (27) 

for m ^ 1. For cr ^ 0, the operator Ce[ph] + OT m Ah is bounded, self-adjoint 
and positive definite and therefore invertible. Since for U m G Hh the right- 
hand side of equation (27) also belongs to Hh, we find that the eqution has a 
unique solution U m G Hh- 



Equation (26) with the help of property (25) implies the first energy equal- 
ity 



M 

I / 1 \ I I i r\ t t m I I ^ II t t/'/t'I m II 9 



m=l 

M 



Ce[p h ],oo ~ r II 1 1 C e [c ft ], oo_ 

A/ 

+ 11^^^)111^ = i||^||^ ]!0O + ^(F-,^)-Kr m , (28) 

m=l 

compare with (17). Similarly to the proof of Proposition 1, it implies bound 
(24). 

Remark 2 Proposition 3 remains valid for any ph, bu and Ch satisfying 
the conditions ph^p>0,bh^v>0,Ch^O and sup (phj + bhj + Chj) < oo. 

m 

Remark 3 The solution U = U a to the finite- difference scheme (21)-(23) 
(specified in Proposition 3) depends continuously on o ^\ for any 9 ^ \. 

Indeed, let o~q ^ 1/2. Then the difference U a — U ao satisfies the operator 
equation in 

Ce[p h ] d t (U a - UJ + A h (U a - U ao )^ = (a - a)A h (U ao - UJ on u 7 . (29) 

Let I h Vj := Ei^j-i V kh+i/2 for j ^ 0. Since -11^ = and A h W = 
d x {—bhd x W + IhCo[ch]W) onuh,oo, the property follows from bound (24) (with 
FW = 0) applied to equation (29). 

We define the mesh counterparts of the norm in L 2 (X, oo) such that 



oo 



h 

W\\l h :=^W} + £ W?h, \\Wf bh := £ Wfh, 

j=J+l j=J+l 

OO 

W\\i,D h ■■= (stW)jWjh+ (*eW)sWfc 

j=j+i 

concerning the correctness of the last definition, see [11]. 

Corollary 1 Let FJ 1 = and U®j = forj ^ J andm ^ 1. If the solution 
jjm ^ j Qr a n m ^ to the scheme (21)- (23) satisfies the approximate 



2 

se,D h 



TBC (10) with some operator S = «S re f, then an equality 

M 

-b^Y,^^l)uf )m T m = \ PoQ \\U 

m=l 

M 

+ E - ^rnPoo\\d t U m \\l Dh + 6ooP^ Wm ||| ft + Cooll^lll^J ^ 
m=l 

holds for all M ^ 1. /is right-hand side is nonnegative for a ^ \. 

Proof. By virtue of equation (21) at the node xj with FJ 1 = 0, relation 
(10) is equivalent to the boundary condition (8); thus, the solution to the 
scheme (21)-(23) satisfies the scheme (6)-(9) as well. Taking the difference 
of the energy equalities (28) and (17) (with G = 0) and applying simple 
identities 

n^iia = ml + \\m h , w\\l h = \W\\l h + liwillv 

\\w\\% M , x = \\w\\% M + x h j\\w\\l tDh 

for any W G H^ and x/j, = ph-,Chi we obtain the announced equality. 
We also derive stability in another norm. 

Proposition 4 Let F = F<® + d x F^ with F^ m G H h and F^ m G H h 

for any m ^ 1 and 11® G H^. Then, for a and < j, the second energy 
bound 

maX {oS/ l|[/1 ^^ l|f/|1 ^} ^ 



iiii 



2 M 



+ \& Fmm ^ + TM^ 1 ' ^ + ^ fW mm hj™) (30) 



c p 

m=l 



holds for the solution U m G Hh, for all m ^ ; to the finite- difference scheme 
(21)-(23) and any AO 1. Here 



\U 



|2 

l(2),oo 



^a-\)r\\d t U\\c Hh6 



+ \\\\d t U 



tv \\Ce\ph],oo 



M 



The bound holds also in the case 6 = \ provided that F^ = 0. 



Proof. The second energy equality 

M M 



?77.= 1 TO=1 

M 



\\\Ul\?c Hh , + Y.( F ' n ' d ^ m )u,J^ < 31 > 



holds, compare with (20). Similarly to the proof of Proposition 2, it implies 
bound (30). 

Corollary 2 Under the hypotheses of Corollary 1, an equality 

M 

-boo £ (££u?) 9 f c/jv m = § (fcocp^ll^ + cooll^li;^) 

m=l 

M 

+ £ [ax>P*£HI;U + - D^m (b^idxdtu™^ + cboii^ii;^)] r m 

m=l 

holds for any M ^ 1. its right-hand side is nonnegative for a ^ ^. 

Proof. The result is derived by taking the difference of (31) and (20) (with 
G = 0). 

By definition, the discrete TBC is an approximate TBC (10) with the oper- 
ator S = <S rc f . It will be explicitly constructed in the next section. Corollaries 
1 and 2 clarify the energy meaning of conditions (14) and (18) for the discrete 
TBC, for $ = Uj, and are exploited below to prove the conditions (for any 

*)• 



5 Derivation and analysis of the discrete TBC 

Now we turn to derivation of the explicit form for the discrete TBC in 
the form (10) and verification of inequalities (14) and (18) for it. We confine 
ourselves by the case of the uniform mesh U r , i.e., r m = r for m ^ 1. Consider 
an auxiliary finite-difference problem on the uniform part of the infinite mesh 



m x 



pooSedtU + A h ,ooU i<T) = on (uj hj00 \uj h ) x lu t , (32) 
U\ j=J = $ : with 1$!^ :=supg - m |$ ro | < oo, <D° = 0, (33) 

[7° = for (34) 
for some qo ^ 1. Here the limiting finite-difference operator 

A h ,ooW := -&oo§r + CooS^W on W ftj00 \ W A . 
has appeared. We seek for the solution satisfying the following property 

H/7112,00,, := supgW V \U™\ 2 ) < oo (35) 

for sumcently large q > qo. 

Since the coefficients are constant and the meshes are uniform, the stated 
problem can be solved explicitly. For a mesh function u T — > C such that 
I ^ |oo, < 00 f° r some g > 0, recall the reproducing function 



TO=0 



i.e., analytic in the disc Di/ q := {\z\ < ^ C C, satisfying a bound 

\$(z)\ ^ W°°> q for |z|<g _1 . (36) 
1 — q\z\ 

Conversely, for a function p G A(D r ) (for some r > 0), the transformation 
<3> = T~ l [p] such that 

p(m)(Q) 1 /* 27r p(^ 



for any m ^ 0, < 7*1 < r, (37) 



ml 2tt Jo z m 
is well defined implying the Cauchy inequality 

l$loo,i/n < max|p(*)|- (38) 

Hereafter z is the imaginary unit, and Rez and Imz are real and imaginary 
parts of z G C. 



Taking into account conditions (34) and (35), for \z\ < r = i, we calculate 

T pooSedfUj + Ah,ooU) (z) = poose — 
L J J r 

= d(z)(t? j+1 (z) - 2 7 (z)C? i (z) + ^-1(2)) (39) 

provided that 7^ 0, with := cr + (1 — a)z. Hereafter, for j ^ J — 1, 
we extend £7™| •= so that <9+C/? = ?7- cr ' )0 = 0. The coefficients j(z) and 
d(z) are expressed by formulas 

= 2ai0(z - 1) + (1 - 2a ^ (fT) , 7 (*) := 1 + — ~ Z) + ^ 



d(z) 



h Poo n h Cqq 

ai := — — > 0, a := -7— ^ 0. 

2T&0O 2&OC 



By virtue of (32) and (39) a difference equation 

U j+1 (z) - 2 7 (z)^(z) + Uj-^z) = for j > J (40) 
holds. The corresponding characteristic equation has the form 

z/ 2 (z) - 27(2)1/(2) + 1 = 0. (41) 
Notice that d(0) = for a = cr := 

Lemma 1 For cr > 0, a ^ cr anrf # ^ | ; £/ie quadratic equation (41) has 
roots v\,vi £ A(D r ), for sufficiently small r > 0, swc/i t/iat 

z/ 1 (z) = Zf 1 ( 7 (z)), 0<M2;)|<1, ^(z)^^))^, M*)I>1, 



where Z^ 1 ^) = 7 + (— l) fc+1 ^/ 7 2 — 1, k = 1,2 are analytic branches of the 
two-valued inverse function to the elementary Zhukovskii function Z(z) = 
^(z + z -1 ) defined in C w?t/t £/ie cross-cut along the segment [—1, 1] of the real 
axis [16]. 

Proof. The presented formulas are rather elementary. The property 
1/1,1/2 £ ^4(A0 holds provided that 7 (-D r ) C C \ [—1,1]. For validity of 



the latter property for sufficiently small r > 0, it is required that |7(0)| < oo 
(i.e., o 7^ do) and |7(0)| > 1. For o 7^ o"o and 2clq6 7^ 1, formulas 

ai + a a (l-4fl)ai + [2 + (l-4fl)a ]<7 

7U (1 - 2a e)(a - aoY 71 j + (1 - 2a 9)(a - a ) 

hold. Since a > and 9 ^ 4, the nominators of the both formulas are 
positive and thus 7(0) > 1 for d(0) = (a - cr )(l - 2a #) > 0, or 7(0) < -1 
for d(0) = (a — c )(l — 2ao#) < 0. If 2ciq6 = 1, then once again 7(0) = 
1 — do — a\o fa\ < —1. 

The following result corresponds to Proposition 5.3 in [11]. 

Proposition 5 For a > 0, a 7^ cr and ^ 7, t/ie solution to the problem 
(32)- (35) exists, is unique and is given by a formula 



ui = r 1 



vt J (z)<S>(z) 



for j > J - 1. (42) 



This solution satisfies a bound 

\\U\\2,oo, q <C|$| 00, g 

(43) 

/or sufficiently large q > go- -^or read «'£ ?s read too. 

Proof. Let a 7^ ctq. By taking into account Lemma 1, for z G D r the 
general solution to the difference equation (40) has the form 

Uj(z) = Cl (zy- J+ \z) + c 2 ( Z y- J+ \z) for 

with any c\(z) and 02(2;). By virtue of bounds (35) and (36) we find that 
c^iz) = 0, and then from condition (33) we derive a formula (taking into 
account that v\{z) 7^ 0) 

U :j ( z ) = $( z )vf J {z) for j^J-1. (44) 

Since Qv{~ J G A(D r ) exploiting Lemma 1, if the solution to the problem 
(32)-(35) exists, then it is given by formula (42). 

Conversely, the function given by formula (42) satisfies an equation 

T pooSedtUj + A,oo^j a) {z) = for j ^ J, z G D r 

and therefore equation (32) as well. It also satisfies conditions (33) and (34). 



By virtue of (38) and (36) we get that, for any j ^ J and m ^ 1, bounds 



gr m |C/ri<maxK- J ( Z )$( Z )KC^^ 

Ul=ri 1 — (/o r l 



hold for gi = j- > - > qo, with 

Co(J-i) = ^lmin = min |z/i(z)|, C j = u lmax = max |^(z)| for j ^ J. 



Therefore bound (43) holds with C = , Cl , where 

°° 2 1 1 

C l = S ^o/ ' = ~2 + 7—2 ' Cl > °" 

j = j_\ lmin lmax 

For real the functions $(2) and 7(2;) are real as well for z G R. If in 
addition \z\ ^ r, then \"f(z)\ > 1 and thus ^i(z) = Z-f (7(2:)) is real. Therefore 
U is real too, see (37). 

The proofs of Lemma 1 and Proposition 5 remain valid also for o = 0, 
< \ and 6 ^ 0. 

We go back to the derivation of the discrete TBC. By virtue of formula 
(44) we have 

U j+1 (z) - U^{z) = L(z) - Uj(z) for j > J- 



r 



dJuW- — (poodtU + CoJjW)) 



±d{z)fa-i*)(z)*(z) (45) 



Therefore it is easy to check that a formula 
6h 2 

/ oo ' J 

holds, compare with (39). By virtue of the well known formula for the multip- 
klication of two poer series it leads to the discrete TBC (10) with the operator 
S = <S re f of the discrete convolution form 

sz^ m := ^ (r * *r = ^ E Rq ® m ~ q for m > ( 46 ) 

q=0 

with the kernel 

R-T- 1 [d(z){v 1 -v 2 )(z)]. (47) 
Let us see that Propositions 1 and 2 on stability are valid for S = <S rc f. 



Proposition 6 For a ^ ^, a 7^ cr and ^ |, £/ie operator S = «S re f 0/ t/ie 
discrete TBC (46) satisfies inequalities (14) and (18). 

Proof. We apply a method first suggested in [9] . Fix any M ^ 1 and real 
values $ x , . . . , $ M . Extend $ m = for m = -1, and m > M. We define a 
function U by formula (42) for j ^ J — 1 and set, for example, UJ 1 := for 
^ j < J - 1 and m ^ 0, and then set F := C^p/J^C/ + A^ (ff) - 

By virtue of Proposition 5, the constructed function U serves as the real 
solution to the problem (32)-(34) and therefore as one to the scheme (21)- 
(23), where F = on (coh,oo\^h) x <^ r and = 0. Then Corollary 1 implies 
inequality (14) whereas Corollary 2 implies inequality (18) since U™ = $ m 
for ^ m ^ M. 

The case a = 00 ^ 2 w ^ n be a ^ so covere d in Remark 4 below. 

We find the kernel R explicitly. We introduce quantities by recurrence 
formulas 

2 7T7. — 1 777/ — 1 
Pm,a,/3 = PPm-l,a,p CtPm-2,a,P for m ^ 1, (48) 

m m 

P0,a,p = 1, Pm,a,P = for 171 < (49) 

with parameters a, /3. Their close connecton to the classical Legendre poly- 
nomials will be shown below. 

Proposition 7 For a > 0, a ^ ctq or a = ctq ^ \, and 6 < 4, £/ie kernel R 
of the operator of the discrete TBC (46) is real and is expressed by an explicit 
formula 

R m = 2aiV5 \Pm,a,p - apm-2,a,p] for m ^ 0, (50) 

with the quantities a, (3 and 5 of the form 

a = a a 1: p = — - — , (51) 

_ 1 _do_ _ 1 _ rf (l-4fl) + rfi 

ao l + ^o' ai " (l + ^o)(l-4^) + ^i' 1 j 

(J := (1 + ad )[(l + <jd )(l - 40) + trdj > 0, (53) 

, a Coo . n , 2 r 
d := — = — t ^ 0, di := — = 4 — — > 0. (54) 
a>\ Poo 0,1 h l 



Proof. Let first a ^ ctq. By Lemma 1 for z G D r we have 

(v l -v 2 )(z) = 2{/ 1 \z)-l. (55) 

One can straightforwardly verify that 

d\z) ( 7 2 (.) - 1) = 
= a\(l -z + d Q z^)[{l - 40) (1 -z) + (di + do(l - A6))z {<7) ] = 

= a\5{aoz — l){a\z — 1) = a\5{az 2 — 2f3z + 1), (56) 

where the coefficients do, d±, oto, a\ and a, (3, 5 are given by formulas (51)-(54). 
Let first also We write down a formula 

az 2 -2fiz+l = (xz) 2 - 2/ixz + 1, 



where x = Jol and a = -7= are real for a > 0, or x = i\/\a\ and a = —i—^= 
v p ^ VI I p ^ 

are purely imaginary for a < 0; herewith x 2 = a and x/i = (5 are always real. 

Inserting the formula into (56) leads to 

d(z) </j 2 (z) - 1 = -a x V5 {/(xz) 2 - 2/jxz + 1 (57) 

at least for sufficiently small z (in accordance with the proof of Lemma 1), 
where j/w is an analytic branch of y/w on C with the cross-cut along the 
negative real half-axis Kew < such that yl = 1- Formulas (55) and (57) 
imply an equality 

a( \( \( \ o /I (^) 2 ~ + 1 

d(z)(z/i - v 2 )(z) = -2 ai V5 

\J\yczy — 2[i>cz + 1 

The following generalized formula for the reproducing function of the Legan- 
dre polynomials holds [9]: 



Xj^Kzf - 2hkz + 1 



m=l 

for any x G C, integer / ^ and sufficiently small z, that easily follows from 
the classical one in the case x = 1, / = [16]. Therefore 

oo 

- = -2a lV ^^ x m [P m _ 2 (/x) - 2/xP TO -i(/i) + P to (m)] z m , 

m=0 



where P m (fi) = for m < 0. Applying the recurrence relation for the Legendre 
polynomials 

TTL — 1 TTL 

fiP m -i(ji) = 2m _ 1 Pm-2(jj) + 2m -l Pm ^ for m ^ °' ( 58 ) 
one can simplify the last formula as follows 

d(z)(y! - v 2 )(z) = 2 ai V6j2 (Pm - Pwr-*)W 

i.e., to derive a formula 

R m = 2 ai V5-^-—{P m -P m - 2 ){ri, m^O. (59) 
2m — 1 



m=0 



Following [6, 13], we introduce modified Legendre polynomials p m ^(z) := 
x m P m (z). From (58) clearly p m , a ,p = Pm,x{p) satisfy recurrence equalities (48) 
and (49) and, in particular, they are real. Therefore formula (50) is proved. 

For a = formula (57) is simplified and takes the form d(z) \J^f 2 {z) — 1 = 
— a\y/~5y/l — 2(3z. Hence one can easily verify that formulas (48), (49) and 
(50) remain valid and even are simplified in the case a = 0. 

Owing to continuous dependence of U on a ^ ^ ( see Remark 3) and R on 
a ^ | (it is clear), one can pass to the limit as a — > cr on the left-hand side 
of the discrete TBC (10) with S = <S rc f and on the right-hand side of equality 
(46). This justifies the validity of formula (46), for R of the form (50), in the 
case a = o"o ^ \ (that is possible only provided that 4ai | 2ao ^ # < <^)- 

Remark 4 Proposition 6 remains valid also in the case a = cr ^ \. To 
see this, it suffices to insert formula (46) into inequalities (14) and (18) for 
S = <S rc f that have been already proved, for o ^ o~o, and pass to the limit as 
a — > do taking into account the continuous dependence of R on a. 

In practical computations, recurrence equalities for R are more convenient 
than formula (50). 

Proposition 8 The kernel R satisfies the recurrence equalities 

Rm = 2m-3 _ x _ m-3 aRm _ 2 ^ 
m m 

R° = -2aiVS, R 1 = 2 ai V5(3. (61) 



Proof. The form of formula (59) differs from the corresponding one in 
[11], Proposition 5.7 only by a constant multiplier. Therefore Proposition 5.8 
in [11] implies a recurrence formula 

m m 
Inserting x/i = (3 and k 2 = a leads to (60). Formulas (61) straightforwardly 

follow from (50) and (48), (49) for m = 1. 

Typical graphs of lg \ R m \ are presented on Figures 1 for Examples 1 and 

2, see Section 6 below (the values of parameters are given on Figures 2 and 

3). 
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Figure 1: Graphs of lg|i? m | in dependence with m, for 9 = 1/12, in Example 1 (left) and 
Example 2 (right) 

One can rather easily extend the above results are to the case of the third 
boundary condition at x = 0, or to the Cauchy problem where equation (1) 
is posed on R. The case a < 1/2 could be also analyzed under suitable 
additional condition between steps r and hj. 

6 Numerical experiments 

Consider the initial-boundary value problem (l)-(3) for the simplest ho- 
mogeneous heat equation where p{x) — 1, b(x) = 1, c(x) = and f(x,t) — 0, 
for ^ t ^ T. In Example 1, we base upon an exact solution 

Ul ( X ,t) = ./-*°I e -(™. ) 2 /[4(^o)] 
V t + to 

with parameters x* > and to > and take the data g(t) = Ui(0 : i) and 
u°(x) = m{x,0). We choose x* = 1.25, to = 0.03125, T = 1 and X = 2.5. 




Note that |u°(a:)| < 3.8 • 1CT 6 for x ^ X. On Fig. 2, we demonstrate the 
numerical solution and its error computed for a = ^, 6 = h = 0.05 and 

1 1500 - 



x1Cf 6 

1, 3-, 








Figure 2: Example 1. The numerical solution (left) and its error (right) for a = |, 9 = 
h = 0.05 and r = -r^r 

1500 

In Table 1, the absolute errors are given in dependence with M (where 
Mr = T) for 6 = and = tL For = 0, they do not practically change for 
M ^ 500 whereas for (9 = the error continues to decrease up to M = 2000 
thus allowing to reach values of 725 times less. Note that actually the value 
Uj^O has been used and the minimal reached absolute error is less than this 
one. 



M 


20 


50 


100 


200 


500 


1000 


2000 


= 


4.39 • 10~ 2 


8.30- 10~ 3 


1.30 ■ 10~ 3 


5.38 ■ 10" 4 


8.63 • 10" 4 


9.15 • 10~ 4 


9.28 ■ 10" 4 


6 = 5 


4.53 ■ 10- 2 


9.50 • 10- 3 


2.20 ■ 10~ 3 


4.89 ■ 10" 4 


7.28 • lO" 5 


1.35 • 10~ 5 


1.28 • 10~ 6 



Table 1: Example 1. The absolute errors in dependence with M for h = 0.05 



Notice that if one sets simply the Neumann boundary condition (i.e., takes 
S = 0) instead of the discrete TBC at x = X, then it is necessary to increase 
X three times to reach the error of the same order of smallness (for the same 
h and r); herewith the maximum absolute error is reached at (x,t) = (X,T) 
(the corresponding graphs are omitted). 

In Example 2, we take the data g(t) = t 2 and u°(x) = 0. The exact 
solution to such problem U2(x,t) = 32£ 2 /4(£) is also known and is calculated 



by applying the recurrence formulas [17] 



Jo(0 = erfc(0 = 




-£erfc(Q, 



4(0 = ^"4-2(0 " -/n-i(0, n = 2, 3, 4, 
2n n 

where £ = Choose T = 1 and X = 1. 

On Fig. 3 we present the numerical solution and its error computed for 
a = |, 9 = h = 0.1 and r = 0.01. In contrast to Example 1, now the 
solution is not close to for x = X. The error is maximal at the node 
(xj,t m ) = (h,r) (but not on the artificial boundary x = X). Notice that 
decreasing of X down to 0.2 for the same mesh steps does not increase the 
error (for u° = and while applying namely the discrete TBC, this is natural 
and clear from above). Moreover, if once again one sets simply the Neumann 
boundary condition instead of the discrete TBC at x = X, then (for the same 
h and r) the absolute error equals 0.15 and is unacceptably large. It decreases 
to the same values as on Fig. 3 only when X increases five times. 






Figure 3: Example 2. The numerical solution (left) and its error (right) for a = |, 9 = y^, 
h = 0.1 and r = 0.01 

On Fig. 4 we give the graphs of errors in the cases 9 = and 9 = |. Their 
forms are different and the maximum absolute errors are about two orders of 
magnitude greater than in the case 9 = j^- 

In addition in Tables 2 and 3 we put the absolute errors in dependence 
with M for h = 0.1 and various 9. For 9 = 0,^ and | the errors do not 
practically change already for M ^ 50; herewith their minimum values for 
9 = and 9 = | are close whereas for 9 = | one is approximately twice larger. 



Figure 4: Example 2. The error of numerical solution in the cases 9 = (left) and 9 = | 
(right) for a = ±, h = 0.1 and r = 0.01 

In contrast, for 6 = ^ the error continues to decrease rather rapidly up to 
M = 500 thus allowing to reach values about four orders of magnitude less 
(though it increases slightly when M grows further). 



M 


5 


10 


20 


50 


100 


200 


9 = 


0.020 


6.989 • 10" 4 


4.239 • 10~ 4 


3.506 • 10" 4 


3.402 • 10" 4 


3.376 • 10~ 4 


e = J2 


0.020 


4.948 • 10" 4 


1.224 • 10~ 4 


1.799 • 10~ 5 


4.700 • 10~ 6 


9.300 ■ 10~ 7 


° = l 


0.019 


4.654 • 10" 4 


2.565 • 10" 4 


3.237- 10~ 4 


3.340 • 10" 4 


3.366 • 10~ 4 


° = i 


0.019 


4.355 • 10" 4 


5.931 ■ 10~ 4 


6.613 • 10" 4 


6.717- 10" 4 


6.743 • 10~ 4 



Table 2: Example 2. The absolute errors in dependence with M for h = 0.1 



M 


300 


400 


500 


600 


650 


6 = i 


2.638 • 10- 7 


6.797- 10- 8 


4.032 ■ 10~ 8 


8.962 • 10~ 8 


1.062 ■ 10" 7 



Table 3: Example 2. The absolute errors in dependence with M ^ 300 for h — 0.1 and 9 = ^ 
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